Electrically driven spin excitation in a ferroelectric magnet DyMnOs 
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Temperature (5-250 K) and magnetic field (0-70 kOe) variations of the low-energy (1-10 meV) 
electrodynamics of spin excitations have been investigated for a complete set of light-polarization 
configurations for a ferroelectric magnet DyMnOa by using terahertz time-domain spectroscopy. 
We identify the pronounced absorption continuum (1-8 meV) with a peak feature around 2 meV, 
which is electric-dipole active only for the light i5-vector along the a-axis. This absorption band 
grows in intensity with lowering temperature from the spin-coUinear paraelectric phase above the 
ferroelectric transition, but is independent of the orientation of spiral spin plane {be or ab), as 
shown on the original Pa (ferroelectric polarization) || c phase as well as the magnetic field induced 
Pa II a phase. The possible origin of this electric-dipole active band is argued in terms of the large 
fluctuations of spins and spin-current. 

PACS numbers: 75.80.-l-q, 76.50.-|-g, 75.40.Gb, 67.57.Lm 



I. INTRODUCTION 



The magnetoelectric (ME) effect, in which the electric 
polarization can be induced by a magnetic field H or 
the magnetization can be induced by an electric field E, 
has been known since the 1950'sii However, the materi- 
als showing the large ME efi'ect are rare. Recently, the 
interest has been revived by the observation of the ferro- 
electric polarization Pg and its 90° flop in the presence of 
a finite field in orthorhombically distorted perovskite 
manganites, iJMnOs, where R represents Tb, Dy, and Gd 
ionsi^'^'^ In DyMnOa, for example, the coUinear spin or- 
der of Mn ions sets in at the Neel temperature Tj^" of 
39 K with the incommensurate (IC) modulation vector 
g^" ~ 0.36 along the 6-axis4 Noticeably, the fcrroelec- 
tricity shows up along the c-axis when the transverse be 
spiral (cycloidal) spin order— evolves below Tc of 19 K. 
The direction of Pg is fiopped to the a-axis when H is 
applied along the &-axisi^ So far, a wide variety of the 
ferroelectric magnets, termed multiferroics, have been 
found in a family of cycloidal-spin magnetsii The fer- 
roelectricity in i^MnOa can be microscopically explained 
by the spin-current model^ or equivalently by the inverse 
Dzyaloshinski-Moriya effect f^i^ as expressed by the rela- 
tion, Pg oc Bij X {Si X Sj), where e^j is the unit vec- 
tor connecting i and j sites. Namely, the cycloidal spin 
structure can induce the ferroelectricity along the direc- 
tion perpendicular to the modulation vector and within 
the spiral spin plane. A direct evidence for this mecha- 



nism was recently demonstrated by controlling the vector 
chirafity {Si x Sj ) of TbMnOs with 

In the magnetically induced ferroelectric phase, the 
strong ME coupling would also produce the intriguing 
low-lying optical excitation. In PMnOs, there are ex- 
perimental (Ref. [lol ) and theoretica l^^i^^ arguments on 
the possibility of the presence of the low-lying spin ex- 
citation, called electromagnon, which is activated elec- 
trically. Far-infrared spectroscopy has provided the use- 
ful information for the electrodynamics of TbMnOar^ 
GdMnOa^iH and Eui_:rY.^Mn03^ii^ revealing that an 
additional light absorption contribution indeed emerges 
in the ferroelectric phase at terahertz (THz) frequencies 
below the optical phonon frequencies of the perovskite 
structure. Among a family of PMnOs, DyMnOs shows 
the largest ME response; the change of the dielectric con- 
stant upon the Pg fiop exceeds 500% and P, reaches 2 
mC/m^i^ii Therefore, DyMnOa offers a good opportu- 
nity to clarify the low-energy electrodynamics coupled 
with the spin excitation. 

Here we report on a complete set of optical spectra 
at THz frequencies of a multifcrroic DyMnOa with vary- 
ing the light-polarization (both the electric field E'^ and 
magnetic field iJ" components of light), temperature T, 
and applied H, to unravel the strongly electric-dipole 
allowed spin excitation. The observed excitation is con- 
firmed to be active mainly for E'^ \\ a and not to be 
correlated to the direction of the spiral spin plane. We 
focus on spectral signatures, being highly sensitive to the 
light-polarization, in a variety of ME phases that can be 



tuned by both T and H. The format of this paper is as 
follows. In Sec. [IT] we describe the detail of the experi- 
mental procedure for measurements of THz time-domain 
spectroscopy. Section Hill is devoted to the optical inves- 
tigations of the low-energy spin dynamics in a variety of 
ME phases of DyMnOs, in which we clarify the unique 
selection rule of the spin excitation driven by or 
by measurements of the light-polarization (Sec. IIII Ap . 
T (Sec. ImB]) . and H (Sec. HlTC]) dependences. On the 
basis of the results of Sec. Illli the origin of the electri- 
cally driven spin excitation observed at THz frequencies 
in DyMnOa is discussed in Sec. IIVI The summary of 
this paper is given in Sec. El We also describe in Ap- 
pendixes the estimation procedure of the complex opti- 
cal constants from the raw data obtained by THz time- 
domain spectroscopy. 

II. METHODS 
A. Sample preparation and characterization 

The large single crystals of DyMnOs were grown by 
a floating zone method. The available ac, a6, and be 
surfaces of the crystals with a typical sample size of 4- 
16 mm^ were cut from the ingot by using a back Laue 
photograph. The respective samples were polished to 
thickness of 150-800 /im. The polishing procedure fre- 
quently induces the damage on the surface for the case 
of the orthorhombically distorted perovskite manganites. 
We confirmed that the optical response was almost un- 
changed before and after polishing procedure. 

The samples cut from the same ingot used in this work 
were characterized by the measurements of the dielectric 
constant at low-frequencies as functions of T and H. The 
fcrroelectricity along the c-axis in zero H and the Pg flop 
from the c to a-axcs in H were clearly observed around 19 
K and 20 kOe at 7 K, respectively. The obtained phase 
diagrams in a plane of T and H are consistent with the 
previously reported one of DyMnOa A 

B. Setup for terahertz time-domain spectroscopy 

We used the THz time-domain spectroscopy in a trans- 
mission geometry to obtain complex refractive index n 
without using Kramers-Kronig (KK) analysisJ^i^i THz 
time-domain spectroscopy can easily access the low- 
energy electrodynamics in solids, typical frequencies from 
0.1 THz to 3 THz.— The lower limit of the measurable 
frequencies depends on the size of the samples since 0.1 
THz corresponds to 3 mm in wavelength. We adopted a 
transmission geometry and the sample was attached to 
the perforated Cu holder with a diameter of 2 mm, 2.8 
mm, 3 mm, and 3.5 mm, which was selected to suitably 
match the size of the samples. 

Figure [H shows the schematic illustration of THz time- 
domain spectroscopy in a transmission geometry in this 
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FIG. 1: (Color online) Schematic illustration of the experi- 
mental setup for THz time-domain spectroscopy in a trans- 
mission geometry. 



study. For the measurements, we used the photoconduct- 
ing (PC) sampling techniquoi^ to obtain THz radiation 
pulse. The femtosecond laser pulses delivered from the 
mode-locked Tiisapphire laser (center wavelength of 800 
nm, pulse width of 90 fs, and repetition rate of 80 MHz) 
were used as a source. The laser pulses were divided by 
a beam splitter (BS). One was used as the pump and 
another was used as the gate. The pump pulses were ir- 
radiated on to the ZnTe crystali^ to induce the optical 
rectification effect, acting as a source of THz radiation.— 
The wire grid (WG) polarizers were inserted in between 
off-axis paraboloidal mirrors to obtain the linear polariza- 
tion. The polarized THz pulse was focused on the sample 
by off- axis paraboloidal mirrors. The cryostat was placed 
within the box filled with dry nitrogen to eliminate the 
absorption of water. The gate pulses were introduced to 
the low-tcmperature-grown GaAs (LT-GaAs) photocon- 
ducting device coupled with a dipole antenna after the 
appropriate time delayi^ The pump pulses were mechan- 
ically chopped and the photocurrent induced by the elec- 
tric field of the THz pulse was lock-in detected. There- 
fore, the induced photocurrent, which depends on both 
amplitude and phase of THz radiation, can be obtained 
by varying the optical delay line. Further details of our 
scheme for THz generation and detection can be found 
in Refs. [l^ andfisl. 

For the measurements of the H dependence, the sam- 
ple was placed in a cryostat equipped with a super- 
conducting magnet. H supplied from a superconduct- 
ing magnet was applied along the 6-axis up to 70 kOe. 
The measurements in H were mainly performed at 7 K. 
We adopted a Faraday geometry. Further details of our 
experimental setup combined with the superconducting 
magnet can be found in Ref. [23 . 
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FIG. 2: (Color online) Low-energy electrodynamics of a complete set of available configurations of single crystals of DyMnOa, 
measured around 10 K in zero magnetic field. Upper and lower panels show the real and imaginary parts of e/x spectrum 
(closed circles) when i5" and ff" were set parallel to the crystallographic axes of ac, ab, and be surface crystal plates. The solid 
lines in (a) and (b) are results of a least-squares fit to reproduce low- and high-lying peak structures by assuming two Lorentz 
oscillators for e. On the contrary, spectra in (d), (e), and (f) can be reproduced by two Lorentz oscillators for e and fi, as 
indicated by solid lines. 



Since we consider the spin excitation, which can be 
driven by either E'^ or H'^, we numerically estimated h 
and used the quantity of n'^ = efJ,{uj) (not = e), where 
e is the complex dielectric constant and fi the complex 
magnetic permeability. Details of the estimation proce- 
dure and the validity of the quantities of e/z to character- 
ize the spin excitation are described in Appendixes. 



III. RESULTS 

A. Light-polarization dependence 

First, in Fig. [2l we show e/z spectra (closed circles) in 
the be spiral spin phase {Pg \\ c) of DyMnOs to determine 
the selection rule with respect to the light-polarization, 
in which E'^ and H'^ were set parallel to the crystallo- 
graphic axes using available ac, ab, and be surfaces of 
the crystal. Upper and lower panels of Fig. [5] present, 
respectively, the real (Re) and imaginary (/m) parts of 
efi spectra, measured around 10 K with H = 0. There 



are remarkable peak structures and noticeable optical 
anisotropy with respect to E'^ and H'^ . Such a light- 
polarization dependence, including the negligible absorp- 
tion for E"^ \\ e and iJ" || 6 [Fig. [^c)], as well as the 
conspicuous spectral change around Tc and Tj^", cannot 
be ascribed to the crystalline-electric-field excitation of 
/ electrons, which was reported for TbMnOa around 4.7 
meV by inelastic neutron scattering studies i^^'^^ 

A broad and pronounced absorption continuum (1- 
8 meV) overlapped by two peak structures around 2 
mcV and 5.5 meV is observed in Im[eii] spectrum for 
E'^ II a and iJ" || c [lower panel of Fig. [2[a)]. In ac- 
cord with this, a dispersive structure can be discerned 
in Re[efi\ spectrum [upper panel of Fig. [UJa)]. Such a 
low-lying peak structure was also identified previously 
at 2.9 meV for TbMnOg^ 2.5 mcV for GdMnOs^i^iH 
and 3 meV for Eui-^rY-cMnOa (although this value de- 
pends slightly on x).— Noticeably, /m[e/i] of DyMnOa 
reaches 10 (corresponding to the optical conductivity of 
about 3 cm~^), which is roughly a factor of 2-5 
larger than the absorption of TbMnOa and GdMnOs, 
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FIG. 3: (Color online) Temperature dependence of the real (upper panels) and imaginary (lower panels) parts of the selected 
e^L spectra (closed circles), of DyMnOa for (a) i?" || a, (b) H'^ \\ a, and (c) H'^ \\ c. Note that the scale of the vertical axis 
is different in respective figures. The Im[efj,] spectra shown in (b) and (c) are vertically offset for clarity. We also plot the 
extracted Im[e] and Im[fi] spectra [at (b) 5 K and (c) 4.2 K] indicated by dotted and gray lines, respectively, by assuming that 
Im[efj,] spectrum can be represented by two Lorentz oscillators for e and jj, (solid lines). The Im[fj] spectrum is multiplied by 
5. On the contrary, efi spectra for || a [(a)] can be fitted by two Lorentz oscillators. 



suggesting DyMnOs as the prototype in the study of the 
low-lying excitation. Such a strong absorption includ- 
ing the spectral shape, peak position, and magnitude, 
was clearly confirmed to be observed only for || a, 
as evidenced by the e/j, spectra for E'^ || a and i/" || b 
[Fig. HUb)], i.e., with different H'^ polarization. The low- 
temperature spectrum can be fitted with two Lorentz os- 
cillators, as indicated by solid lines. The ratio of Re[efj] 
for i?" II a to other configurations was about 2, being 
consistent with the dielectric constant Re[e] at 10 kHz4i^ 
This indicates the large anisotropy extending up to THz 
frequencies, associated with the electric-dipole active ex- 
citation. 

In Im[en] spectrum for E'^ \\ c and H"^ \\ a [Fig. HJd)], 
on the other hand, a broad peak feature is observed where 
the peak energy nearly matches the low-lying peak en- 
ergy for i?" II a. However, the magnitude of /m[e/i] in 
this configuration 0.3) is much smaller, by a factor of 
30, than the magnitude of Im[efi] for E'^ \\ a. By com- 
paring the spectra for E'^ \\ c and iJ" || a [Fig. [D^d)] 
and E"^ \\ b and H"^ \\ a [Fig. [DJe)], this absorption ap- 
pears to arise from the H'^ \\ a stimulation, in accord with 
the conventional magnon selection rule for the be spiral 
spin statCfi^ We also discern the broad peak feature in e/i 
spectrum for E'^ \\ b and H'^ \\ c [Fig. ^f)] whose magni- 



tude of /m[e/i] 0.3) is comparable to those of /m[e/i] 
(~ 0.3 — 0.5) for cases of \\ a stimulation. We cannot 
distinguish the selection rule of this tiny absorption on 
the basis of the light-polarization dependence alone due 
to the presence of the intense absorption for i?" || a and 
H"^ II c [Fig. ^a)]. As will be shown in Sec. Hirei how- 
ever, we identify that this sharp peak structure is also 
driven by H'^. 



B. Temperature dependence 

We present in Fig. ^h) the T dependence of the se- 
lected e/i spectra for 77" || a. With increasing T, /m[e/x] 
tends to decrease in magnitude and the peak position is 
shifted to the lower energy from 2.8 meV at 5 K to 2 
mcV at 36 K [the lower panel of Fig. Mh)]. Above T^", 
the broad peak structure is diminished, as can be seen in 
Im[e^] spectrum at 50 K (a finite value of Im[eij] comes 
from the absorption tail of the optical phonons). Ac- 
cordingly, there is no dispersive structure around 2 meV 
in Re[efi] spectrum [upper panel of Fig. [3jb)]. We also 
show in Fig. (SJc) the selected e/i spectra for || c as 
a function of T. With increasing T, the peak position is 
also shifted to the lower energy from 3.3 meV at 4.2 K 
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FIG. 4: (Color online) Temperature (T) dependence of the 
integrated spectral weight per Mn-site, A^efi , as defined by Eq. 
U]) in the text, of the ac surface crystal plate of DyMnOa for 
II a [(a)] and \\ a [(b)]. The upper panel of (a) shows 
the phase diagram with variation of T. PE and PM stand for 
paraelectric and paramagnetic, respectively. We also show the 
schematic illustrations of the be spiral and collinear spin states 
induced by T. Neel temperature T^"^ (39 K) and ferroelectric 
transition temperature Tc (19 K) are indicated by vertical 
solid lines in (a) and (b). The data measured at 75 K and 
249 K were also plotted at the 50 K position in (b) with square 
and triangle symbols, respectively. The horizontal dashed line 
in (a) represents the estimated contribution of the optical 
phonon absorption at 5 K. 



to 2.3 meV at 31 K and the magnitude of /m[e/i] also 
tends to decrease. The observed T variation is similar to 
the case for iJ" || a; the peak structure seems to van- 
ish above T^" = 39 K. Even though we cannot identify 
the clear signature of the peak structure above 31 K in 
/m[e/x] spectrum [the lower panel of Fig. [2{c)] due to 
the detection limit in our experiments 1.8 meV), as 
determined by the size of the sample, there is no signa- 
ture of the dispersive structure around 2 meV in Re[e^ 
spectrum at 46 K [the upper panel of Fig. [Sljc)]. On 
the basis of the T dependence as well as of the light- 
polarization dependence presented in Sec. IIII Al it is rea- 
sonable to consider that the observed peak structure can 
be assigned to the spin excitation of Mn ions driven by 
H'^ II a or TJ'^ || c. Indeed, in inelastic neutron scattering 
experiments on the be spiral spin phase of TbMnOs , the 
magnon dispersion along the [100] or [001] axis comes 
across fc = around 2 meVj ^'^'^^ To estimate the con- 
tribution of the magnon to /m[e/i] spectrum, we plot in 
Figs. [3^b) and[3Jc) the Iml/i] spectra for H'^ || a at 5 



K and for H'^ \\ c at 4.2 K, respectively, as indicated by 
solid lines. These were derived from the relation, /m[e/x] 
(= Re[e\Im[ii\ + /m[e]i?e[/x]) by assuming Lorentz oscil- 
lators of e and ^ for contributions of phonon (a dotted 
line) and magnon, respectively. 

Next, we show in Fig. El^a) the T dependence of the 
e/^ spectra for E'^ \\ a. At 249 K, only the absorption 
tail due to the optical phonon is observed. With decreas- 
ing T but even above T^", it is clearly discernible that 
Im[e^] in the low-energy region below 8 meV increases 
and forms a Debye-like absorption band [the lower panel 
of Fig. [3)Ja)]. In accord with this, the dispersive struc- 
ture emerges in i?e[e/x] spectrum [the upper panel of Fig. 
|3i;a)]. In the coUinear phase below T^" 39 K), a 
broad peak structure emerges around 2 meV and grows 
in intensity continuously with further decreasing T. In 
the be spiral spin phase below Tc (=19 K), the additional 
peak structure shows up around 6 meV and its peak po- 
sition shifts to the lower energy while the low-lying peak 
structure shifts to higher energy. 

To be more quantitative, we deduced the integrated 
spectral weight per Mn-site, A^eff, as given by 



N, 



eff 



2moV 



uj'Im[e{uj')^i{Lj')]dLj', (1) 



where toq is the free electron mass, e the elementary 
charge, and V the formula (DyMnOa) cell volume. The 
integrated energy range was chosen between ui = 0.8 
meV and uj2 = 9.5 meV to nearly cover the full band. 
Figure (Hb) presents T variation of the estimated A'eff. 
We also plotted A^off at 75 K (square) and 245 K (tri- 
angle) at the 50 K position for clarity. A"eff for i?" || a 
gradually increases below T]^" = 39 K and is sharply 
enhanced below Tc = 19 K (indicated by vertical solid 
lines). This is in contrast to A"off for H'^ \\ a (we adopted 
uji = 1.5 meV and uj2 = 4.6 meV to estimate A'off of the 
magnon band), as shown in Fig. [H^a); A'off reaches the 
maximum at Tc where the be spiral spin structure evolves 
from the collinear spin state [note that A^off for H'^ \\ a 
includes the contribution of the phonon structure, as in- 
dicated by a dashed line in Fig. Ufa), which could be 
crudely estimated with the /TO[e] spectrum at 5 K, as in- 
dicated by a dotted line in Fig. El^b)]. The electric-dipole 
active {E'^ \\ a) excitation almost disappears above Tj^" 
[Fig. [31Ja)] and its energy range (0-10 meV) nearly coin- 
cides with that of the magnon band, which was revealed 
by inelastic neutron scattering experiments on the be spi- 
ral spin phase of TbMnOa i^i^^ This coincidence of the 
spectral energy region also ensures that the observed ex- 
citation is magnetic in origin. 



C. Magnetic field dependence 

As an origin of the electromagnon, the rotational mode 
of the spiral spin plane around the 6-axis (propagation 
vector) has been theoretically proposed based on the 
spin-current modelf^ This was taken into consideration 
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FIG. 5: (Color online) 
Magnetic field (H) depen- 
dence of the real (up- 
per panels) and imaginary 
(lower panels) parts of the 
selected e/x spectra (closed 
circles) of the ac surface 
crystal plate of DyMnOa 
at 7 K when (a) E'^ or (b) 
H"^ was set parallel to the 
a-axis. The ferroelectric 
polarization flop occurs at 
20 kOe ior H \\ b. Note 
that the Im[e^] spectra in 
(b) are vertically offset for 
clarity and that the scale 
of the vertical axis is dif- 
ferent in (a) and (b). 



in the analysis of the inelastic neutron scattering re- 
sults for the be spiral spin phase of TbMnOa^ and the 
observed dispersion relation (2.4 mcV at fc = 0) was 
assigned to the rotational mode of the be spiral spin 
plane. These apparently contradict the results of e spec- 
tra of GdMnOa pSiii^ in which the peak structure with 
Im[e] ~ 1, comparable to Im[e] of TbMnOa, was also 
identified in the IC spin-coUinear paraclcctric phase. In 
the course of the Pg flop from the c to a axis at He ^ 
20 kOe, IC g^" (~ 0.37) is kept nearly constant for 
DyMnOsr^ eis contrasted by the IC to commensurate 
transition for TbMnOa.— Therefore, the measurements 
of the effect of H on e/i spectra of DyMnOs can provide 
further insights into the nature of the electromagnon. 

The upper and lower panels of Fig. El^a) present the 
selected i?e[e/i] and /m[e^] spectra for i?" \\ a, measured 
at 7 K for varying H with H \\ 6-axis, respectively. The 
experimental setup with use of a superconducting mag- 
net (up to 70 kOc) restricted the measured energy range, 
and thus we concentrated on the behavior of e/i spectrum 
below 5 meV. Even though the low-lying peak structure 
grows in intensity and its peak position is shifted to the 
lower energy above He, the broad continuum-like absorp- 
tion scarcely depends on H {He « 20 kOe at 7 K); the 
i?" II a polarized absorption is similarly identified as the 
most notable electric-dipolc active band even in the ab 
spiral spin phase. Such a tendency is also seen in H de- 
pendence of TVoff, as shown in Fig. ^h). We selected 
Ui = 0.7 meV and UJ2 = 5 meV to estimate A^cft of the 
low-lying excitation. The negligible effect of H is clearly 
discerned even H exceeds He- This behavior is different 
from the reported result of TbMnOa, in which the low- 
lying electromagnon disappears by application of H up 
to 80 kOe at 12 K (Ref. In the experiment of Ref. 
[lol ). H was applied along the c-axis, thus giving rise to 



the suppression of Ps along the c-axis, according to the 
ME phase diagram of TbMnOa,^ Namely, the flop of 
from along the c-axis to the a-axis or equivalently the 
rotation of the spiral spin plane from be to ab cannot be 
induced unless the H is applied along the a or 6-axis, as 
we adopted here. In addition, we measured the e/i spectra 
for E'^ II a at 17 K just below Tc in H along the 6-axis up 
to 70 kOe, to confirm the same tendency. Furthermore, 
we gathered the e/t spectra for E'^ || c and H'^ \\ a, at 
7 K in up to 70 kOe. The upper and lower panels of 
Fig. [S^b) shows the selected Re[efi] and Im[efi] spectra, 
respectively. The estimated iVoff is also presented in Fig. 
[S^a) (we used uji = 2.4 mcV and W2 = 4 meV). Above 
He, the direction of E^ (|| c) becomes perpendicular to 
the spiral spin plane and hence would activate the elec- 
tromagnon arising from the rotational mode of the ab 
spiral spin plane. Even though the evolution of the tiny 
absorption is visible around 2 meV above He, the strong 
absorption, to be compared with /m[e/i] for E'^ \\ a in 
zero H, is not discerned at all for E'^ \\ e, but only for 
E"^ II a. 



IV. DISCUSSION 

In the spiral spin phase of i?Mn03 , the electromagnon 
or equivalently the rotational mode of the spiral spin 
plane is expected to become active when E'^ is set per- 
pendicular to the spiral spin plane, i.e., E'^ \\ a and 
i?" II c in 6c and ab spiral spin phases, respectively. This 
scenario is theoretically considered on the basis of the 
spin-current mechanism (Ref. [T^ that is now recognized 
to explain the magnetically driven ferroelectricity in a 
variety of multiferroics found so far, including RMjiO^^ 
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FIG. 6: (Color online) Magnetic field (H) dependence of the 
integrated spectral weight per Mn-site, A'cff , as defined by Eq. 
11]) in the text, of the ac surface crystal plate of DyMnOa for 
H"^ II a [(a)] and E'^ \\ a [(b)], measured at 7 K. Upper panel 
of (a) shows the phase diagram with variation of H. We also 
show the schematic illustrations of the be and ab spiral spin 
planes induced by H. H of 20 kOe at which the ferroelectric 
polarization flop occurs at 7 K, is indicated by the vertical 
solid line in (a) and (b). The difference of A^cff at 7 K with 
[in (a) and (b)] and without H [in Fig. [4], arises from the 
difference of the integrated energy range (see text). 



Indeed, a single peak structure around 2.5 mcV has been 
identified in the e spectrum for i?" || a in the be spiral 
spin phase of TbMnOs (Ref. [13). Accordingly, one of 
low-energy branches of the magnon band, i.e., around 2 
meV at A; = 0, in the he spiral spin phase of TbMnOa was 
assigned to be responsible for the electromagnon model 
by inelastic neutron scattering experiments (Ref. [25l ). 

In the present study, we could uncover new features 
of the low-energy spin excitation driven by either or 
H'^ in a variety of spin phases of DyMnOa on the basis 
of the detailed light-polarization, T , and H dependences, 
as presented in Sec. IIIII First, our result can provide the 
compelling evidence that the possibility of the rotational 
mode of the spiral spin plane is excluded as the origin of 
the intense absorption around 2-3 meV for E'^ \\ a, that is 
also observed in the be spiral spin phase of TbMnOa we 
firmly identify the electric-dipole active spin excitation 
along the a-axis well survives even in the ab spiral spin 
phase induced by applying H up to 70 kOe along the b- 
axis [Fig. [5ja)] . This observation clearly contradicts the 
former theoretical^ and experimental^^ arguments de- 
scribed above. Second, we clarify that the electric-dipole 
active spin excitation along the a-axis can be regarded 



as the continuum-like absorption band with two broad 
peak structures (1-8 meV) [Figs. [DJa) and^b)] rather 
than the single peak structure (1-5 meV), as reported in 
the be spiral spin phase of TbMnOs (Ref. [13). Third, 
we identify the spin excitation with a single peak struc- 
ture when i/" is set parallel to the a-axis or the c-axis 
[Figs. [Ud)-12I^f )] ; which can be assigned to the antifer- 
romagnetic resonance (AFMR) of Mn spins at fc = 0. 
The k = magnon mode for if^ || a can be reasonably 
ascribed to one of the lower branches of the magnon, as 
observed around 2 meV at fc = in the inelastic neu- 
tron scattering spectrum of the be spiral spin phase of 
TbMnOa, even though it was previously assigned to the 
rotational mode of the spiral spin planci^ This complex- 
ity partly comes form the fact that the peak position 
of AFMR for i?" || a nearly matches the peak position 
of the lower- lying peak structure for E'^ \\ a, as we re- 
vealed here [Figs. Ufa) and[2)Jd)]. On the contrary, the 
observed broad continuum-like band for E'^ \\ a spreads 
over the wide energy range (1-10 meV), which cannot be 
explained by fc = magnon absorption alone. 

The observed broad peak and continuum-like (0-10 
meV) shape of the electric-dipole active band does sug- 
gest that the spin fluctuation related to the excitation 
such as the two-magnon is more plausible rather than the 
Goldstone mode like electromagnon. The optically active 
two-magnon excitations are possible over the full magnon 
bandwidth with the combination of +k and — fc magnons, 
and hence should distribute the spectral weight up to the 
energy twice of the magnon bandwidth. According to 
the inelastic neutron scattering study on TbMnOs (Refs. 
[I^IH), the magnon band width is about 8 meV. This 
value is considerably narrowed as compared with that 
of LaMnOs (- 32 meV)^ or PrMnOg (- 20 meV)^ 
with the A-typc antiferromagnetic ground state, due to 
the larger orthorhombic lattice distortion as well as to 
the strong spin frustration effect in the cycloidal magnet 
TbMnOa. The magnon bandwidth of DyMnOs, in which 
the inelastic neutron scattering study would be difficult 
to perform due to the presence of the neutron-absorbing 
Dy ions, should be comparable to, or slightly smaller 
than, that of TbMnOs, since DyMnOs is structurally 
close to TbMnOs and shows the similar be plane cycloidal 
order at the ground statc4i Therefore, two-magnon en- 
ergy range is anticipated to be up to 10-15 mcV, which 
well fits the observed energy range of the electric-dipole 
active spin excitations. 

Here we propose that the observed unique optical 
dipolc activity may share the same origin, i.e., the spin- 
current mcchanisnif^ with the static Pg. The strong 
E'^ II a selection rule irrespective of the be or ab spi- 
ral spin plane does suggest the large mutual fluctua- 
tion within the ae plane, of the neighboring spins lin- 
ing up along the c-axis, but not along the b nor a- 
axis; namely {6 Si x 6Sj) \\ b, Cij \\ e, and hence SP'^ 
[cx Cij X {5Si X SSj)] II a. The antiferromagnetic ex- 
change interaction along the c-axis in the perovskite man- 
ganite is likely the origin of the large spin fluctuation. 
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FIG. 7: (Color online) (a) Sample geometry with the complex Fresnel transmission coefBcients from vacuum to the sample 
tv-ts and from the sample to vacuum ts^v The complex transmission coefBcient t was obtained by Eq. ([Bill , d and n are the 
thickness and the complex refractive index of the sample, respectively, (b) Typical example of the measured THz wave form in 
time-domain with and without the sample (ac surface crystal plate of DyMnOa), measured at 5 K in zero magnetic field, i?" 
and were set parallel to the c and a-axes, respectively. The amplitude of the transmitted THz wave form with the sample 
was vertically offset for clarity, (c) Power transmission spectrum of DyMnOa. 



contrary to the ferromagnetic coupling along the a-axis. 
This scenario, though a more elaborate theory is needed, 
may qualitatively explain the strong optical activity for 
E'^ II a also in the spin-coUinear paraelectric phase as ob- 
served. Then, broad peaks around 2 meV and 5.5 meV 
may be assigned to the density-of-state singularities of 
the two-magnon band. 



V. SUMMARY 

In summary, we have measured the spectra of a 
complete set of possible configurations of single crystals 
of DyMnOs by using THz time-domain spectroscopy. A 
rich variety of optical excitation at THz frequencies was 
identified by detailed measurements of light-polarization, 
T, and H dependences. In particular, the spin excitation 
driven by E'^ was confirmed to become strongly electric- 
dipole active along the a-axis, irrespective of the be or 
ab spiral spin plane. The broad (1-10 meV) band fea- 
ture indicates the contribution of two-magnon like ex- 
citations arising from the large fluctuations of spin and 
spin-current inherent to this spiral magnet. 
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APPENDIX A: RAW DATA OF TERAHERTZ 
RADIATION IN TIME-DOMAIN 

A unique advantage for the use of THz time-domain 
spectroscopy is to directly monitor the electric field of 
the radiated THz pulse in time-domain.— This tech- 
nique yields the amplitude and phase spectra by the fast 
Fourier transformation (FFT) of the measured THz wave 
forms with and without the sample. Therefore, we need 
no Kramers-Kronig transformation, which is necessary 
for obtaining the complex optical spectra by conventional 
optical spectroscopy. As a typical example, we show in 
Fig. Wljo) the transmitted THz wave form with the sam- 
ple in time-domain, performed at 5.2 K in zero H. The 
sample was the ac surface crystal plate of DyMnOa . E'^ 
and H'^ were set parallel to the c and a-axes, respectively. 
We also plotted the reference of THz pulse used in this 
setup. The amplitude of the transmitted THz wave form 
is reduced due to the absorption of the sample, as com- 
pared to the reference. Its phase is also delayed due to 
the refractive index of the sample. The ratio between 
FFT spectrum with and without the sample, 

%^ = VTexpi-icj)) (Al) 

can yield the power transmission T and phase (f> spectra. 
The example of the obtained T spectrum is shown in 
Fig- [7{c). We can identify the clear absorption peak 
structure around 3 meV. This peak can be assigned to 
the spin excitation driven by iJ" || a, as presented in Fig. 

m 



APPENDIX B: RELATIONSHIP BETWEEN 
TRANSMISSION COEFFICIENT AND OPTICAL 
CONSTANTS 

We numerically derived the complex optical spectra 
without the KK transformation by adopting the follow- 
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FIG. 8: (Color online) (a) Validity of using Eqs. ()B2a|) and 
HB2b|) . instead of Eqs. (|B3a|l and (|B3b|) . The /m[e/i] spec- 
tra of DyMnOs for E'^ \\ c and |t a at 5 K was almost 
unchanged even if we used Eqs. (|B3aTl and (|B3b)) . The fixed 
parameters of Re[^] and Im[fi] are described in the right-hand 
side of the figure, which were confirmed to be reasonable by 
assuming that e can be represented by a single Lorentz os- 
cillator [see, (b)]. (b) The estimated Im[n] spectrum (closed 
circles) by assuming that e can be represented by a single 
Lorentz oscillator. We also plot the Im[n] spectrum (a solid 
line), as extracted from Eq. (|B1|I combined with Eqs. (|B2a[) 
and ((B2bl) . 



ing procedure. Without taking into account the effect of 
the multiple reflections inside the sample, the complex 
transmission coefficient t for bulk in normal incidence 
as schematically shown in Fig. [Tj^a), can be simply ex- 
pressed as 



t = iv— s^s— V cxp (^—i—d{n — l)j = 



E, 



sample 



rcf 



(Bl) 



where lu is the frequency, c the velocity of light, d the 
thickness of the sample, and n (= u + ik) the complex re- 
fractive index of the samplci^ Experimentally, the effect 
of the multiple reflections can be avoided by restricting 
the time range of the FFT. The complex Fresnel trans- 
mission coefficients^ from vacuum to the sample iv->s 
and from the sample to vacuum ig^v, can be, respec- 
tively, given by 



n + 1' 

2h 
n + 1' 



(B2a) 
(B2b) 



when the complex magnetic permeability is 1. In this 
case, n can be expressed as ^/e, where e is the complex di- 
electric constant of the sample. Accordingly, h and thus 
e can be numerically calculated from Eq. (|Bip using ex- 
perimentally determined values; T and (p in Eq. (jAl[) . 
This procedure is widely recognized to estimate h and e 
of a variety of materials including dielectrics, semicon- 
ductors, and metalsj ^^i^^ 

On the contrary, in magnetic materials, n should be 
exactly expressed as ^/ejl near the magnetic resonance, 
which usually locates around ~ meV (Ref. IsTl ). In this 
case, Eqs. (|B2a[) and (IB2bp should be, respectively, re- 
placed by 



2fi 

n + 
2n 
n + n 



(B3a) 
(B3b) 



APPENDIX C: ESTIMATION PROCEDURE OF 
COMPLEX OPTICAL SPECTRA 

In the cases of the measurements for i?" \\ a using 
ac and ab surface crystal plates of DyMnOa, we exper- 
imentally identified from light-polarization dependences 
[Figs. IHa) and [UJb)] that the observed large absorp- 
tion originates from E'^ response. The H'^ response in 
these configurations can be considered to be negligible 
compared to E'^ response. This is confirmed from the 
results using be surface crystal plate where there is no 
absorption in H"^ \\ b (S" || c) [Fig. djc)], and only a 
shght absorption in H'^ \\ c {E'^ || b) [Fig. [^f)]. There- 
fore, we assume /i = 1 in the Fresnel coefficient and used 
Eqs. (|B2aP and (jB2bp . while keeping n ~ y/ejl in Eq. 
(jBl|) for generality. This is also applied to the case for 
E'^ II c and H'^ \\ 5, in which the negligible absorption 
is clearly identified [Fig. HJ^c)]. We numerically solved 
Eq. ((BT|) combined with Eqs. ((B2a| and (|B2bp to match 
the measured values given by Eq. (|Aip and determined 
the e/i spectra. For the sake of convenience, we used the 
quantities of Re[efi\ and Im[efi\ in this paper, although 
they can be precisely expressed as Re[e\ and /?T!,[e] in the 
above condition. 

On the other hand, we identify a peak structure around 
3 meV, which only appears when H'^ was set parallel to 
the a-axis [Figs. [2jd) and[2Ue)] or the c-axis [Fig. [U^f)]. 
These small absorption can be assigned to response. 
In addition, a broad background absorption is observed 
above 4 mcV as seen in the power transmission spectrum 
at 5.2 K [Fig. [7jb)], which can be attributed to the tail of 
the optical phonon lying around 16 meVJ^ In such a case, 
both e and fi could contribute to the transmission spec- 
trum. Accordingly, four independent quantities, i?e[e], 
/m[e], i?e[^], and Im[fj], cannot be obtained only from 
the transmission measurements of T and 0. However, 
in the case of a thick sample, the main signature of the 
transmission spectrum is dominated by bulk absorption 
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and the effect of /i in the Fresnel coefficient is considered 
to be small, as discussed below. Therefore, we used the 
Fresnel coefficient of Eqs. (jB2aP and (jB2b[) . while the ef- 
fect of i?" response was taken into account by adopting 
n = in Eq. (jBip . The validity of this approach was 
firmly confirmed following calculation. 

We numerically solved Eq. (jBip by putting Eqs. (jB3aP 
and (|B3bp with i?e[^] and Im[fi] as fixed parameters with 
respect to oj and estimated Im[efi] spectra of DyMnOa 
for E'^ II c and H'^ \\ a at 5.2 K. For examples, we show in 
Fig. mja) the calculated Im[efx] spectra at selected values 
of Re[^] and Im[^]. The fixed parameters used here are 
provided in the right-hand side of the figure. The Im[efj.] 
spectrum, as derived from Eqs. (|B2aP and (|B2bp [i?e[^] = 
1 and Im[ij] = in Eqs. (|B3aP and (jB3bp ]. is also plotted 
by closed circles with a solid line. As can be clearly seen 
in the figure, no significant difference is observed even 
though /i is included to tv^s and tg^v 

To further check the validity of the used values of Re[fi] 
and Im[^], we assumed that e comes from the contribu- 
tion of the optical phonon structure, which can be repre- 
sented by a single Lorentz oscillator. This approximation 
has been conventionally used to estimate the /x spectra 
in antiferromagnets, as exemplified by the cases of the 



orthofcrritesiai We used the following parameters; the 
resonance frequency, the damping rate, and the dielec- 
tric constant for the high frequency above the measured 
range, were set to 9 meV, 0.5 meV, and 16, respectively. 
With this approximation, the Im[fj.] spectrum was nu- 
merically calculated using Eq. (jBip combined with Eqs. 
(|B3aP and (|B3bp . The result is presented in Fig. Mh) by 
closed circles. We also show the Im[fj] spectrum (a solid 
line), as derived from the Im[efi] spectrum that was esti- 
mated by using Eqs. (|B2a|l and (jB2bp . We assumed that 
the Im[efi\ spectrum (= i?e[e]/m[/i] + /m[e]i?e[/i]) can 
be expressed by two Lorentz oscillators for e and fj,. The 
both data coincide with each other below 4 meV, ensuring 
the validity of using Eqs. (jB2aP and (|B2bp . Noticeably, 
the maximum values of Re[^] and Im[ij] are estimated to 
be about 1.01 and 0.02, respectively, indicating that the 
fixed parameters shown in Fig. [Sja) are reasonable. How- 
ever, there is a discrepancy above 4 meV, which mainly 
arises from the lack of the spectral information above 6 
meV and thus from the inappropriate parameters used 
here to reproduce the optical phonon structure. In this 
work, we numerically solved using Eq. ()Bip with Eqs. 
(jB2a[) and (jB2bp and estimated h (= ^/ejl) spectra. 
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